A Klebsiella pneumoniae NDM-1+ bacteriophage: Adaptive polyvalence and disruption of heterogenous biofilms

Bacteriophage KL-2146 is a lytic virus isolated to infect Klebsiella pneumoniae BAA2146, a pathogen carrying the broad range antibiotic resistance gene New Delhi metallo-betalactamase-1 (NDM-1). Upon complete characterization, the virus is shown to belong to the Drexlerviridae family and is a member of the Webervirus genus located within the (formerly) T1-like cluster of phages. Its double-stranded (dsDNA) genome is 47,844 bp long and is predicted to have 74 protein-coding sequences (CDS). After challenging a variety of K. pneumoniae strains with phage KL-2146, grown on the NDM-1 positive strain BAA-2146, polyvalence was shown for a single antibiotic-sensitive strain, K. pneumoniae 13,883, with a very low initial infection efficiency in liquid culture. However, after one or more cycles of infection in K. pneumoniae 13,883, nearly 100% infection efficiency was achieved, while infection efficiency toward its original host, K. pneumoniae BAA-2146, was decreased. This change in host specificity is reversible upon re-infection of the NDM-1 positive strain (BAA-2146) using phages grown on the NDM-1 negative strain (13883). In biofilm infectivity experiments, the polyvalent nature of KL-2146 was demonstrated with the killing of both the multidrug-resistant K. pneumoniae BAA-2146 and drug-sensitive 13,883 in a multi-strain biofilm. The ability to infect an alternate, antibiotic-sensitive strain makes KL-2146 a useful model for studying phages infecting the NDM-1+ strain, K. pneumoniae BAA-2146. GRAPHICAL ABSTRACT


Introduction
Bacteriophages, viruses that infect bacteria, are presently of significant interest to researchers concerned with the proliferation of antibiotic-resistant genes in the environment and the increasing need for alternative therapies (Domingo-Calap and Delgado-Martínez, 2018). This interest is justified, in part, due to the aggressively antimicrobial nature of bacteriophages, which are responsible for 20% of the bacterial turnover on Earth (Suttle, 2007). Bacteriophages destroy antibiotic-resihreferstant bacteria with the same efficacy as any other host and are self-propagating, increasing their numbers exponentially after a low initial dose. Bacteriophage therapy has been proven effective in treating life-threatening infections where traditional antibiotics have failed (Huang et al., 2022) with phage-based products commercially available that include brand names such as PhageGuard (for the control of Listeria) or Intralytix (for treatment of E. coli infections). Despite such examples and increasing support from researchers and doctors who recognize the potential of phage therapy, obstacles to its widespread implementation persist (Sulakvelidze et al., 2001;Loc-Carrilloand and Abedon, 2011;Brüssow, 2019). Limitations primarily include (1) difficulty in getting phage to the relevant host sites at reasonable efficiency, (2) the large numbers of virions that may be required to have an effect, and (3) a generally limited virus host range (Sulakvelidze et al., 2001;Duckworth and Gulig, 2002).
The last obstacle, bacteriophage host specificity, is also considered one of the advantages of phage therapy in that it is possible to target specific problem bacteria while sparing other, perhaps helpful, bacterial species. However, limiting phage infectibility to a single species also becomes a challenge when one needs to apply phage therapy to multiple pathogenic bacterial species or multiple bacteria co-present in environmental applications (e.g., against biofilms or filamentous bulking). Even within the same species, bacterial variations can preclude the efficacy of phage treatment, necessitating precise diagnosis along with large existing phage banks (Nilsson, 2014). For applications where a broad host range is advantageous, such as when confronted with diverse community structures present in biofilms, the search for polyvalent bacteriophages (phages that can infect hosts of different genera or species) is on the rise.
Although biofilms can form from a single species of bacteria, they are most often composed of a diverse range of organisms that include different species of bacteria, algae, fungi, and protozoa (Wimpenny et al., 2000;Balcázar et al., 2015). In complex heterogeneous communities such as this, biofilms are most effective in providing a physical barrier for microbes against environmental stressors like antibiotics and phage predation (Hall and Mah, 2017). Cells within biofilms also tend to be less metabolically active, rendering many antimicrobials ineffective (Stewart, 2015). These are just a few reasons bacterial infections and engineering challenges resulting from biofilm formation are so difficult to address. With antibiotics and disinfectants proving limited in their ability to clear biofilms, bacteriophage present a promising therapeutic alternative.
This manuscript presents a lytic bacteriophage, KL-2146, isolated from a wastewater treatment plant (WWTP) in Salt Lake City, Utah. Its host, Klebsiella pneumoniae BAA-2146 (Hudson et al., 2014), was originally isolated from a US hospital patient in 2010 who had recently received medical care at a clinic in India. K. pneumoniae BAA-2146 carries multiple antibiotic resistance genes (ARG) on multiple plasmids, including the broad range ARG New Delhi Metallo-β-lactamase 1 (NDM-1), making it resistant to most known antibiotics (Yong et al., 2009). Phage KL-2146 was determined to be polyvalent in that it can infect another Klebsiella host species, K. pneumoniae 13,883, which is antibiotic-sensitive. It was also demonstrated to effectively disrupt biofilms possessing multiple Klebsiella strains in contrast to another isolated non-polyvalent phage species capable of infecting only a single strain.
Detailed studies of potentially useful bacteriophages benefit from experimentation with gene knockouts, substitutions, fusions, tags, and targeting endonucleases, e.g., CRISPR (Chen et al., 2019). Because these methods primarily involve the use of antibiotics as selective markers in the host strain, conducting these experiments in bacterial strains that are not susceptible to antibiotics is very difficult. For this reason, and because of the inherent risks in handling such bacteria, the ability of a phage such as KL-2146 to infect multiple hosts (polyvalence), including antibiotic-sensitive hosts, is a valuable feature. This report provides the first example of polyvalence in a phage infecting an NDM-1 positive Klebsiella bacterial strain. It also discusses an observed switching in host infection efficiency between the host strains K. pneumoniae BAA-2146 and K. pneumoniae 13,883.

Host characteristics and growth procedures
Klebsiella pneumoniae ATCC# BAA-2146, possessing NDM-1 (Hudson et al., 2014), is a Gram-negative, non-motile, encapsulated, rod-shaped bacterium. This strain is resistant to most known antibiotics (Yong et al., 2009), with multiple narrow and broad-range antibiotic resistance genes present on several large plasmids. The second potential host that KL-2146 is tested against is K. pneumoniae subsp. pneumoniae ATCC# 13883 (Daligault et al., 2014). This strain possesses the same basic structural features as the NDM-1+ strain, is cultured using the same methods, and has no known antibiotic resistances. Other strains used for determining the polyvalence of phage KL-2146 include K. pneumoniae ATCC# 35657 (Ayu et al., 2014), Klebsiella oxytocca CCUG# 15788 (Singh et al., 2016), which is a close relative of K. pneumoniae, and clinical isolates of K. pneumoniae (19C5A and 19C6A) obtained from Dr. Matt Mulvey's lab at the University of Utah (Mulvey and Simor, 2009). Growth curves of K. pneumoniae BAA-2146 and K. pneumoniae 13883 were determined experimentally from single colony isolations grown overnight in LB (Miller, 1972) at 37°C. Newly 1:100 inoculated cultures were shaken (190 RPM) at 37°C with 1 mL of the sample taken for spectrophotometric measurements (OD 595 ) every 10 min. Colony forming units (CFU/mL) at exponential growth was determined by chilling samples on Ice for 10 min, serial dilution in LB, and plating on LB + 0.7% agar. Detailed growth curves were not established for the other potential hosts, which did not yield plaques from phages used in this study.

Phage isolation, propagation, and physical characterization
Phage KL-2146 was isolated from activated sludge obtained from a wastewater treatment plant (Salt Lake City, Utah) by infection of K. pneumoniae BAA-2146 utilizing methods in sequential filtration and centrifugation previously described (Delfan et al., 2012;Bhattacharjee et al., 2017;Hockett and Baltrus, 2017). A single plaque was thusly isolated and propagated by infection of the same hosts in liquid culture (Earnshaw et al., 1976). Phage lysates from cultures were shaken with chloroform and centrifuged for 20 min at 6000 g. The cleared lysate was then titered for downstream infections of larger (500 mL) cultured using an MOI = 10. Upon lysis, these cultures were concentrated by overnight centrifugation (15,302 g, JA 16.25 Beckman rotor) and resuspended in SMG buffer (50 mM Tris-HCl pH = 7.5, 100 mM NaCl, 10 mM MgSO 4 , 0.01% gelatin). Purification of phages was performed by CsCl isopycnic centrifugation (Earnshaw et al., 1976) using protocols described in Gilcrease et al. (2020). Phage concentrate was layered onto a step gradient of SMG + 10% Sucrose, SMG + 1.4 g/mL CsCl, and SMG + 1.6 g/ mL CsCl in an SW41 Ultraclear™ centrifuge tube for SW41 Beckman rotors. These preparations were then centrifuged at 220,000 g for 2 h at 20°C. The band containing phage, located at the 1.4/1.6 CsCl interface, was extracted with a syringe and dialyzed against SMG buffer using a Spectra/ Por™ 7 50,000 MWCO dialysis membrane. DNA was extracted for sequencing using the Norgen Phage DNA Isolation kit (cat# 46800). Proteins of mature virions were visualized with 12% acrylamide SDS-PAGE gels using the Mini-PROTEAN II gel system, a protein sizing standard (Thermofisher Benchmark Protein Ladder. Cat#10748010), and procedures described in Eppler et al. (1991). Basic phage structural morphology was determined by negative staining transmission electron microscopy (TEM, model JEOL JEM-1400) using methods described in Gilcrease et al. (2020). A phage lacking polyvalence was also isolated to compare phage infection efficiencies toward a multi-strain biofilm. It was isolated on K. pneumoniae 13,883 using the previously described methods and only makes plaques on that strain. This non-polyvalent strain, called phage KL-13883, was not further characterized in this study. Detailed one-step growth curve analysis for phage KL-2146 [n = 3] was conducted by infection of 100 mL host culture at exponential growth phase (OD 595 = 0.5) at 37°C, with a multiplicity of infection (MOI) of 10 phages per colony forming unit (CFU). This higher-thantypical MOI ensures that all bacteria are infected (Nale et al., 2021) and is the same MOI used for large-batch phage purification. This produces total lysis after the first round of infection and allows ease of determining the Plaque Forming Units (PFU)/cell using the CFU/mL known at the time of infection from previous bacterial growth analysis. Two 1 mL samples were then taken every 10 min, one of which was vortexed with 100 μL chloroform. Both were then centrifuged for 5 min at 6000 g, transferred to fresh Eppendorf tubes, serially diluted, and tittered using the agar overlay method (Hockett and Baltrus, 2017). Growth curves were plotted by dividing the PFU yield by the bacteria concentration in culture, which was determined experimentally to be 1.2 × 10 8 CFU/mL at OD 595 = 0.5 ( Figure 1) to obtain PFU/cell. In other experiments where infection efficiency is compared, an MOI =0.026 was used, and samples were taken after only 70 min of infection for comparative analysis. This is when detailed growth curves suggest maximum lysis has occurred after the 1st round of replication in culture. Plaque forming units (PFU) were visualized and quantified using the agar overlay method.

Phage KL-2146 genome sequencing, assembly, and annotation
After extraction of viral DNA, quality was determined by 16S rDNA PCR amplification (to ensure no host DNA contamination) and quantified by spectrophotometry (Thermo Scientific. NanoDrop 1,000). An Illumina TruSeq DNA Library preparation (450 mean insert size) with UDI was performed, followed by Illumina MiSeq 150-cycle paired-end sequencing at the University of Utah HighThroughput Genomics Core Facility. Approximately 29, 303 paired-end reads were trimmed based on a quality score (Phred score) of 28 and de novo assembled using software (Geneious™). GeneMarkS (Besemer et al., 2001) predicted genes were used to search the NCBI non-redundant (nr), conserved domain, COG, and InterPro databases and were annotated using Blast2GO 2.5.0 (Conesa et al., 2005) with accompanying manual curation. ARAGORN (Laslett and Canback, 2004) and tRNAScanSE (Schattner et al., 2005) were used for detecting tRNA genes. Phage KL-2146 mod (KL-2146 re-isolated on the 13,883 strain) was prepared, de novo assembled, and analyzed using the same procedures after trimming 33,212 paired-end reads. Comparative analysis of Frontiers in Microbiology 04 frontiersin.org sequence data for phages KL-2146 and KL-2146 mod was done using Snapgene™. Taxonomy was informed using the Basic Local Alignment Search Tool (BLAST) of the 959 bp major head gene (H1O04_gp48) against NIH sequence databases.

Switching of host specificity and examining polyvalence
To investigate polyvalence and the observed switching of host specificity, phages were separately isolated from KL-2146 infections of different Klebsiella hosts. Phage KL-2146 that was propagated on K. pneumoniae BAA-2146 from individual plaques is referred to as KL-2146. KL-2146 mod was made by propagating KL-2146 through K. pneumoniae 13,883. In making KL-2146 mod , K. pneumoniae 13,883 was used for the isolation of an individual plaque of KL-2146 (MOI = 0.26). This plaque was allowed to diffuse into 1 mL LB overnight and then used to infect 50 mL of K. pneumoniae 13,883 to produce a titered stock. This stock was later used for making larger, cesium-purified phage preparations in the 13,883 strain. In subsequent experiments, the purified stock of KL-2146 mod , which infects K. pneumoniae 13,883 with nearly 100% infection efficiency, was used to infect its original host, K. pneumoniae BAA-2146. Similar to the KL-2146 mod isolation, plaques from the resulting infection were picked and used to reinfect the original K. pneumoniae BAA-2146 culture. The new phage stock is called KL-2146 new . The PFU/mL values given for these experiments are averages taken from [n ≥ 3] experimental repeats.

Infection of multi-strain biofilm
Heterogeneous biofilms containing the two primary K. pneumoniae strains were grown using two different methods. In the first method, biofilms were grown on glass microscope slides using methods adapted from Bhattacharjee et al. (2015). To compensate for the slower growth rate inherent to 13,883, more of this strain was used to inoculate the cultures than the faster-growing BAA-2146 strain. Hence, 25 mL of LB was inoculated with 1:100 overnight cultures of K. pneumoniae 13,883 and 1:250 K. pneumoniae BAA-2146 in a 50 mL conical tube. 75 mm × 25 mm glass slides were inserted into culture tubes and incubated in the upright position overnight at 37°C. Media was replaced by carefully transferring the slides to a tube containing fresh LB every 24 h for 4 days.
At this time, bacteriophages KL-2146, KL-2146 mod , KL-2146 new , and KL-13883 were added to 25 mL of fresh LB to a final concentration of 1×10 7 in four separate 50 mL conical tubes. A fifth tube was prepared without bacteriophages as a control. Glass slides containing biofilm were submerged in the media + phage for 12 h at 37°C. Slides were then washed by dipping them carefully in PBS solution. For CFU quantification of the two strains present in biofilm after 6 days, just prior to staining, samples of biofilm were removed for CFU analysis by scraping across the center of the slide using a 1 cm wide metal spatula. The biofilm substance on the spatula was then transferred to 1 mL of ice-cold LB in a 1.5 mL Eppendorf tube. The samples were serially diluted and plated onto two agar+ LB plates, one with 50 μg/mL ampicillin and the other with "no-drug. " Plates were incubated overnight at 37°C, and colonies were counted the next day. Immediately after the above sampling of biofilm from the slides, slides were dipped into a methanol bath for fixation, air dried, and stained using a Gram staining kit (Brown and Brenn Cat# KTBBR) following the instructions provided therein. Photographs of the biofilm were taken using a light microscope at 100X magnification.
For the second method, biofilms were grown using a custom-made drip flow reactor. Glass microscope slides were seeded for 6 h at room temperature with 25 mL of a 1:100 host culture in 50 mL conical tubes. Slides were then placed at a slight angle into a drip flow reactor so that LB media was allowed to drip over slides at a rate of (1.5 mL/h) for 4 days using a peristaltic pump. Slides were then washed in PBS and stained using a Baclight Bacterial Viability kit (Invitrogen Cat# L7012) following instructions contained in the manual. Biofilms were then visualized using a Nikon Widefield microscope at 60X magnification.

Transmission electron microscopy
For TEM visualization, 1×10 11 PFU/mL of CsCl purified phages were dialyzed in SMG buffer. 3 μL of this mixture was laid onto Formvar 200 Mesh Copper grids that had been subjected previously to a 20 s plasma ionization treatment. After 1 min, the excess phage solution was removed with filter paper. 3 μL of 1% uranyl acetate was added to the grid, left for 45 s, and extracted with filter paper. All negatively stained TEM images were obtained using a TEM model JEOL JEM-1400.

Statistical analysis
The averages of at least 3 technical repeats were taken for CFU and PFU values given, with standard deviation bars or error percentages provided. To determine statistical differences in the number of plaques in growth curve results, one-way analysis (ANOVA, Excel) was performed between the two different data sets. Significance was determined at an alpha level of 0.05 (p ≤ 0.05) and is the standard set for all [n ≥ 3] data that is presented.

Structural and growth characteristics of isolated lytic phages and their hosts
In all experiments, the bacteria strain K. pneumoniae 13883 was observed to grow slower than BAA-2146, lagging behind BAA-2146 in reaching exponential phase by ~30 minutes, as shown in the growth cycle in Figure 1A. K. pneumoniae 35657, K. oxytocca, and sequenced clinical isolates K. pneumoniae 19C5A and 19C6A had growth rates similar to K. pneumoniae BAA-2146 (data not shown).
All of the KL-2146 isolations used in this study (KL-2146, KL-2146 mod , and KL-2146 new ) share identical virus morphology as determined by negative stain TEM; hence only images of KL-2146 are provided ( Figure 1B1). In TEM, KL-2146 was structurally observed to be a member of the Drexlerviridae family (formerly Siphoviridae) possessing an isometric head approximately 60 nm in diameter and a ~150 nm long, filamentous, non-contractile tail. SDS-PAGE analysis of purified KL-2146 virions ( Figure 1B2) showed seven observable structural proteins. KL-2146 and KL-2146 mod (the same phage isolated on different strains of K. pneumoniae) do not show any differences in the size or the total number of protein bands in SDS-PAGE. Phage KL-2146 is shown in one-step growth curves to have a burst size of approximately 142 pfu/cell with a generation time of about 30 min, as shown in the phage growth cycle in Figure 1C. It is a virulent (or lytic) phage based on clear plaque formation and a lack of genes found in sequence analysis that might suggest a lysogenic lifestyle present in the assembled sequence, such as integrases or ParA/B genes (Emerson et al., 2012).

Phage genome sequencing and analysis
To characterize the genetic composition of KL-2146 and to conduct phylogenetic analysis, whole genome sequencing of the purified phage was conducted via Illumina, with the complete genome deposited in GenBank under the accession number MN379832. The KL-2146 genome is 47,844 bps long and 50.8% G + C with circular permutation. Seventy-four encoded proteins are annotated in the genome, including a recognizable large terminase, major capsid protein, several tail proteins including tail tape measure, lysis proteins, and nucleic acid metabolism proteins including DNA helicase, primase, single-strand DNA binding protein, exonuclease, and DNA methylase. Genome sequence analysis revealed phage KL-2146 to be a member of the Webervirus genus within what had previously been described as the T1-like cluster of tailed phages (Casjens and Grose, 2016) or Tunavirinae.

Phage infectivity to its parent host and other strains
KL-2146 was shown to be highly infectious to its host BAA-2146 (the original host from which it was isolated), with complete lysis of culture visible after 30 min of MOI = 10 infection or 90 min of MOI = 0.026 infection. To assess polyvalent properties, phage KL-2146 infectivity was tested using a variety of K. pneumoniae strains phylogenetically similar to BAA-2146 (Supplementary Table 1). Using the plaque overlay method, a faint clearing of bacteria was noted in spot tests on K. pneumoniae 35,657 at low dilutions, but no plaques were produced. No faint clearing from concentrated lysates or plaques was observed upon attempts to infect strains of K. pneumoniae clinical isolates 19C5A, 19C6A, or the K. oxytocca strain. Of the strains tested, only K. pneumoniae strain 13,883 showed clear evidence of infection (plaque formation). KL-2146 is thusly considered a polyvalent phage in that it can infect different strains, or serovars, K. pneumonia. However, plaque titers of KL-2146 on K. pneumoniae 13,883 produced significantly fewer PFUs (approximately 1.6%) than when titered on the original host K. pneumoniae BAA-2146. Seventy-minute infections of BAA-2146 and 13,883 were then conducted in parallel using the original cesium chloride density gradient purified phage KL-2146 that had been isolated using BAA-2146 as a host. In these liquid culture infections, those phages propagated on BAA-2146 were shown when infecting 13,883 to produce an even lower yield (0.17%) compared with infections of BAA-2146 at t = 70 min (Figure 2). The difference in observed yield between direct plaque assay of phage and Frontiers in Microbiology 06 frontiersin.org plaque assays conducted after infection of liquid culture remains unclear, but may be due to (1) more favorable exposure conditions (e.g., more time for attachment) under conditions of direct titration or (2) differences in initial burst size would not be observable on direct titration since lower (or stunted) phage yields from the first round of infection will produce the same number of plaques overnight as higher/normal yields. To assess whether or not KL-2146 infection efficiency would remain stunted after isolation from 13,883 or become acclimated to the new host and infect normally, plaques taken from the less efficient 13,883 infections were further propagated on 13,883 in liquid culture. After concentration and purification, the resulting phage (now called KL-2146 mod ) was found to have acquired normal infection efficiency against this new 13,883 host strain, with PFU yields similar to KL-2146 infections of BAA-2146 ( Figure 2). The KL-2146 mod (with 13,883 as its new host) phage that had been thusly "acclimated" to host 13,883 was shown later to have reduced specificity towards its original host, BAA-2146. This was observed in experiments where phage KL-2146 mod, again purified from infections of 13,883, was seen to produce low phage yields (~15.4%) when infecting the original BAA-2146 strain (Figure 2). In order to determine whether the specificity of KL-2146 mod towards the BAA-2146 strain can be re-acquired after having been lost, the lysate from phage KL-2146 mod was titered on the original host BAA-2146. Plaques were isolated and further propagated on the BAA-2146 host, with the resulting stock now called phage KL-2146 new . This was then used to infect BAA-2146. After 70 min, PFU from infections of BAA-2146 with the KL-2146 new lysate indeed shows a reversal to total infectability of the BAA-2146 strain (Figure 2).

Infections of a multi-strain biofilm
A heterogeneous biofilm composed of bacteria from two separate strains of K. pneumoniae was constructed to demonstrate the comparative efficiencies of polyvalent virus infection versus non-polyvalent infection when two host strains are present. Biofilms were developed on microscope slides using both immersion and drip flow methods and visualized using a Gram Stain and Baclight Viability Kit, respectively. Biofilm heterogeneity was confirmed by counting CFU from the uninfected control on ampicillin plates (1 × 10 7 /mL ± 3×10 6 ) for BAA-2146 and drug-free plates (1.8 × 10 7 ± 2×10 6 ) for the total number of bacteria. This suggests an approximate 1:1 ratio of BAA-2146 versus the antibiotic-sensitive strain 13,883, given the margins of error. In biofilm infection experiments, relative to the uninfected biofilm control, slides exposed to polyvalent phages were mainly cleared of biofilm due to infection of both strains (Figure 3). Backlight Live/Dead color labeling suggests that the comparatively few bacteria remaining on the slides after infection using the polyvalent phage are dead. There was no observable difference between polyvalent isolates KL-2146, KL-2146 mod , and KL-2146 new regarding their efficacy in clearing the multi-strain biofilm. The biofilm slide exposed to the non-polyvalent phage, KL-13883, Six-day heterogenous biofilms of K. pneumoniae 13,883 and BAA-2146 were infected with 1×10 7 phages/mL for 12 h. Top row: Slides seeded with bacteria and submerged fresh media (replaced daily) in were gram stained and visualized at 100X. Bottom row: 60X magnification of biofilm from drip flow reactor (same order as above) stained with BacLite. Green = live, Red = Dead. Biofilms disruption demonstrated on all slides except the control and after infection with the non-polyvalent phage KL-13883. CFU values [n = 3] are from approximately 6 cm 2 sampled from slide surfaces resuspended in 1 mL LB. Error percentages are the standard deviation/CFU. showed substantially less clearing than the uninfected control due to its ability only to infect the K. pneumoniae 13,883 strain. Surviving CFU values were nearly identical to uninfected control, implying a rapid replacement of 13,883 with the resistant bacterial strain. This demonstrates that the polyvalent strain is much more efficient at clearing multi-strain biofilms than a non-polyvalent strain, regardless of the less efficient initial infection of one host or the other.

Phage structure and genome
At the time of this analysis, there were 133 complete genomes of Weberviruses in the public database, the majority of which fit within one of a dozen subgroups of the family previously called Tunavirinae (Korf et al. (2019) and our unpublished results). Figure 4 shows examples of seven of these subgroups, separated by bold horizontal lines. KL-2146 fits in the subgroup called "subcluster D" (Grose and Casjens, 2014) or the "Webervirus" genus (formerly the KP36like genus), all but one of which infects K. pneumoniae (Mishra et al., 2012;Hoyles et al., 2015). Among KL-2146's closest relatives are bacteriophages Sanco (MK618657), MezzoGao (MF612072), GML-KpCol1 (MG552615), KPN N141 (MF415412), KL-2146PN1 (KT001920), KpKT21phi (MK278861), F20 (NC_043469.1) and NJR15 (MH633487), which were isolated in Texas, Maryland, Turkey, Korea, Ireland, Israel, Russia, and China, respectively (Nir-Paz et al., 2019;Richardson et al., 2019;Zurabov and Zhilenkov, 2019;Horváth et al., 2020). This closely related group of Klebsiella phages, all of which are >95% identical in nucleotide sequence, appears to be common and exceptionally widespread across the Earth. The KL-2146 genome reveals a variety of features that make it especially adaptable to host defenses. KL-2146, like many of the Webervirus phages, shares a notable lack of restriction sites such as HindIII, NsiI, NcoI, PaeI, and SnaBI (Hockett and Baltrus, 2017) and expresses its own DNA adenosine methyltransferase (KL-2146_2) suggesting a similar restriction evasion strategy to that of coliphage T1 (Auer and Schweiger, 1984). In addition, a well-studied klebsiella phage, phage KP36, encodes a tail-spike that has a capsule-targeting depolymerase (AXI63_gp50) that is capable of degrading the bacterial capsule and biofilms (Majkowska-Skrobek et al., 2016). The KL-2146 tail protein H1O04_gp67 is 90.32% identical to that protein.

Phage infection efficiency toward hosts
The adaptable host specificity observed, whereby optimal efficiency of infection seemed to switch back and forth between hosts, is interesting and warranted further exploration. Bacteriophages evade host defenses and acquire new host specificities using a variety of strategies. These include the evolution of receptor-binding proteins in tail-spike or tailfiber regions that allow recognition of new hosts, DNA methylation, and DNA sequence modification that protects against host restriction Dot matrix analysis of phage genome DNA alignments within the T1-Cluster of phages, constructed with a word length setting of 13 using Gephard (Krumsiek et al., 2007). The sequences, taken from GenBank, show the diversity of phages within the T1-cluster (Previously known an Tunavirinae). The continuity of diagonal lines indicates sequence similarity between phage genomes. The host genera are color coded. Subgroup D contains the nearest relatives of KL-2146 in the Webervirus genera.
Frontiers in Microbiology 08 frontiersin.org enzymes and CRISPR-acquired immunity (Labrie et al., 2010). Several attempts were thusly made to analyze the mechanism for this reversible switching in host infection efficiency from experimental data.
(1) Restriction Avoidance: Restriction modification (RM) systems in bacteria are known to interfere with phage infection. In these systems, bacterial DNA is modified by methylation of specific sequences that prevent host endonuclease action on that sequence while allowing the targeted degradation of unmethylated foreign (e.g., viral) DNA that may enter the cell, thereby aborting infection (Gregory et al., 2010). In these cases, a phage may attach efficiently to a host but, lacking the host-derived DNA methylation pattern, will be unable to successfully complete infection. This can be overcome at low frequency when a rare infection event allows incoming viral DNA to escape DNA restriction and acquire the methylation pattern after a single round of infection, thereby producing progeny that are insensitive to host RM systems (Labrie et al., 2010). This is consistent with the observations of KL-2146 in Figure 2 regarding its back-and-forth host specificity and is the most likely explanation. The diversity of phages, restriction-modification systems, and host cells makes it difficult to provide a general frequency of restriction avoidance. Studies of this phenomenon using restriction-sensitive phage Lambda mutants suggest this typically happens at a lower frequency, 0.1-0.00001% (Arber, 1965;Rambach and Tiollais, 1974) than observed here for KL-2146 (1.6%). The robust KL-2146 restriction avoidance mechanism (i.e., lack of restriction sites in the genome) may increase the frequency of this avoidance beyond expectations, however.
(2) Differences in Host Growth Rate: While in all infectivity experiments, the bacterial strain K. pneumoniae 13,883 was seen to grow slower than BAA-2146 ( Figure 1A) substitution that is not predicted to affect host specificity ( Figure 5). In summary, the mechanism for host switching appears to be most similar to restriction avoidance, given its behavioral consistencies with what is known about the phenomenon and considering that other mechanistic candidates have been ruled out.

Conclusion
Phage KL-2146 is shown capable of infecting both strains BAA-2146 and 13,883 with equal efficiency only after a single step of inefficient propagation in the other strain, making it a potentially useful model system for the study of phage infection of an NDM-1+ pathogen. The study would have benefitted from including more polyvalent phages infecting NDM-1 for comparative study. It was observed, however, in isolation of these phages from wastewater that there were substantially fewer phages that infected the strain NDM-1 than the antibioticsensitive strain and even fewer that infected multiple strains, including the NDM-1 strain. Indeed, KL-2146 was the only successfully isolated polyvalent phage infector of K. pneumoniae NDM-1+. This is perhaps informative as to the abundance of this type of antibiotic-resistance bacteria present in the wastewater isolation source (Yahya et al., 2015).
The proliferation of K. pneumoniae NDM-1+ strains and other broad-range antibiotic-resistant bacteria in selective reservoirs (such as hospitals and wastewater treatment plants) is an immediate concern, with multidrug resistance genes regularly making it into WWTP effluents and the larger ecosystem (Alexander et al., 2020). Bacteriophages are potentially an effective and ecologically safe antimicrobial tool, one of the few remaining capable of controlling such populations of bacteria in WWTPs. In addition, they can do so in a targeted manner, sparing those bacterial populations necessary for treatment processes. More study is needed to effectively exploit the antimicrobial nature of phages for actual use in such applications. Viral polyvalence is a property that not only allows infections of hosts outside of a single bacterial strain and enhanced removal of biofilms but also allows the detailed study and genetic manipulations of phages infecting these bacteria, particularly those with broad-range antibiotic resistances. This would otherwise be very difficult due to the lack of available selective markers (e.g., antibiotic resistance genes) usable in the resistant host. By exploiting polyvalence and conducting research efforts in an antibiotic-sensitive strain, modern DNA editing methods may be used to construct bacteriophages with enhanced infectivity, drug delivery capabilities, and DNA editing mechanisms such as CRISPR that can target and eliminate specific host genes (Gholizadeh et al., 2020). The study and characterization of phages like KL-2146 are thus an important addition to the emerging field of environmental phage therapy.

Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/genbank/, MN379832.

Author contributions
AB conducted the initial isolation and sequencing of KL-2146. AB and EG conducted propagation/purification of phages and growth curve studies. SC, AB, and EG conducted phage genomic analysis. EG designed the polyvalence study and performed sequence and DNA analysis of KL-2146mod and structural analysis of both phages. EG and AB drafted the manuscript, which RG reviewed. All authors contributed to the article and approved the submitted version.